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Modelling SCUBA sources in a ACDM cosmology: hot 
starbursts or cold extended galactic dust? 
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ABSTRACT 

Previous modelling has demonstrated that it is difficult to reproduce the SCUBA 
source counts within the framework; of standard hierarchical structure formation mod- 
els if the sources are assumed to be the high-redshift counterparts of local ultra- 
luminous infrared galaxies with dust temperatures in the range 40 - 60 K. Here, we 
show that the counts are more easily reproduced in a model in which the bulk of the 
sub-millimetre emission comes from extended, cool (20 - 25 K) dust in objects with 
star formation rates of 50 - 100 Moyr"^. The low temperatures imply typical sizes 
of ~ l(S'85o/lniJy)^/^ arcsec, a factor two to three larger than those predicted using 
starburst-like spectral energy distributions. Low dust temperatures also imply a ratio 
of optical/UV to 850-/im flux which is 30 - 100 times smaller, for the same optical 
depth, than expected for objects with a hot, starburst-like SED. This may help explain 
the small overlap between SCUBA sources and Lyman-break galaxies. 

Key words: submillimetre - galaxies: high-redshift - galaxies: starburst - galaxies: 
evolution - galaxies: luminosity function - infrared: galaxies 



1 INTRODUCTION 



Studies of galaxy properties in the infrared (IR) have added 



their stars (e.g. iRowan- Robinson et al 


I1997I: iMeurer et al. 


19971: iDickinsonI Il998l: iLutz et al.l 


199d: iPettini et al 


19981: iFlores et al.lll999l:lMeurer. Heckman & Calzettlll999l 


Steidel et alJ 'igQgY At low redshifts. 


the bright end of 



the infrared galaxy luminosity function is dominated by 
strongly starbursting systems, which show a rapid in- 
crease in space density towards larger redshifts (see re- 
view by Sanders & Mirabel 199^. The large number of 
bright submillimetre (sub-mm) sources that were detected 
by the SCUBA instrument on the JCMT nevertheless came 
as a surprise ||Smajlj.Jvison_fc Blai 1 I1997I: iHuehes et alJ 
| l998l: iBa raer. Cowic fc Sanders' '19991: iBlain et alJ Il999bf 
Eales et al. 2000; Scott et al. 2002). If the the sub-mm 
sources are the high-redshift equivalent of the strong star- 
bursts observed at low redshifts, the high number counts 
imply that strong evolution in the space density of bright 
IR sources would have to continue beyond 2; ~ 3 - 4. 

Many authors have devised phenomenological models to 
reproduce the sub-mm source counts. These models assume 
a local luminosity function of the sources and a parametrized 
law for the evolution of the sources in luminosity and/or 



density that is constrained to fit the counts at various wave- 
lengths (see lRowan-Robinsonll200ll and references therein). 

It is nevertheless important to understand the ob- 
served sub-mm counts within the standard cosmologi- 
cal paradigm. Previous modelling based on semi-anal ytic 
jGuiderdoni et alJll998l: iDeyriendt fc Guiderdoiiil I2OOOII or 
gas dynamical (ilWda^talll20^l) simulations, has encoun- 
tered considerable difficulties in reproducing the brightest 
sub-mm sources. If the bright SCUBA sources are the high- 
redshift counterparts of the local ultra-luminous infrared 
galaxies (ULIRGs), their dust temperature should be high 
(40 - 60 K). They are then inferred to be forming stars at 
rates of a few hundred to a thousand solar masses per year, 
the exact value depending on the detailed assumptions re- 
garding extinction, initial mass function and the duration 
of the burst. Such extreme star formation rates are very 
difficult to achieve in standard hiera r chical galaxy forma- 
tion models. IPevriendt fc GuiderdonJ (|20o3) show that un- 
der these assumptions, they are only able to reproduce the 
sub-mm counts with a model in which essentially all baryons 
present in massive dark matter haloes at redshifts ~ 3 turn 
into stars over timescales of order 10* years. They propose 
that high merging rates in the early Universe may be re- 
sponsible for this extremely efficient conversion of gas into 
stars, but d o not attempt to m odel this in detail. On the 
other hand. Ifilain et al.l lll999al) have modelled the merg- 
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ing rates of dark matter haloes in hierarchical cosmologies 
and find that, in order to fit the counts, the halo mass-to- 
infrared light ratio of a typical galaxy merger must be 200 
times smaller at redshift 3 than at the present day. This is 
cleaxly an extreme requirement. 

The form of the spectral energy distribution (SED) of 
the sub-mm galaxies is, however, a s ource of major un - 
certainty in the predictions. Recently iFardal et alJ (1200 ll) 
modelled the sub-mm counts using N-body plus smoothed 
particle hydrodynamics simulations and showed that bet- 
ter agreement with the data was obtained if the dust tem- 
peratures were at the lower end of the range usually as- 
sumed for starbursting gala xies. This follows suggestions by 
|^owai^Robinso^j200l|) and lEfstathiou fc Rowan- Robinsm] 
i2002. hereafter ERR^ that the sub-mm emission from high- 
redshift galaxies could arise from cirrus-like emission of 
cold extended dust. An extended distribution of cold dust 
may, for example, be caused by wind-driven dust outfiows 
llHeckman. Armus fc MilevI I199(TI : lAlton. Davies fc Bianchil 
Il999^ . 

In this paper, we explore the degeneracy between the 
inferred star formation rate and dust temperature explicitly 
and treat temperature as a free parameter when fitting the 
sub-mm counts. We demonstrate that hierarchical galaxy 
formation models can reproduce the sub-mm counts with 
moderate star formation rates if the sub-mm emission arises 
from cool, extended dust and discuss the implications for 
the size of the emit ting region. We use the galaxy formation 
mo del described inlKauffmann. White fc Guiderdonl lll993ll 
and iKaufFmann et alTi^QsIl and the diffuse cirrus emission 
SEDs of ERR. 

In 121 describe the galaxy formation model used to 
predict the star formation rates and outline the difficulties of 
hierarchical models in accommodating a high space density 
of strong starbursts out to large redshifts. discusses our 
modelling of the SED and explores the degeneracies in the 
dust models. In we present our results and a comparison 
with observations and in i|S]we discuss their implications for 
the nature of the sub-mm sources. presents our conclu- 
sions. 



2 THE GALAXY FORMATION MODEL 

2.1 Star formation in hierarchical galaxy 
formation models 

We have used the star formation ra tes predicted by the sem i- 
analytic galaxy formation model of lKauffmann et al.l il993l) . 
Here, we briefiy recapitulate some of the assumptions of the 
model relevant for our predictions of the sub-mm counts. 
The model and semi-analytic prescriptions are described in 
detai l in the above reference and elsewhere l|Lacev fc Cole , 
| l993l: iKa uffmann et all Il999l: ISomerviUe fc Primackj Il999l: 
ICole et aL .2000.1 . The calculations are for a fiat ACDM 
model with f^A = 0.7 and Ho = 70 kms"^Mpc"^ The 
model follows the formation and evolution of galaxies within 
a merging hierarchy of dark matter haloes. Simple prescrip- 
tions take into account the dissipative cooling of the gas, the 
formation of stars, feedback from supernovae and the merg- 
ing rates of galaxies located within the same dark matter 
halo. Star formation is modelled by the following equation: 



z = - 

z = 0.5 

0-1 Fi__, z = 1.2 - 




M / Mg yr"' 

Figure 1. The combined star formation rate function at various 
redshifts. 

M = a- , (1) 

where M is the star formation rate, A/com is the mass of 
the available cold gas, tdyn is the dynamical time of the 
galaxy, and a is a parameter controlling the efficiency of 
star formation. Mcoid depends on the coohng rate of the hot 
gas in the dark matter halo, the star formation rate and the 
efficiency of supernovae in heating the gas. In the model, a 
is a free parameter that is set to reproduce the luminosity 
of a Milky Way-type galaxy. The star formation efficiency is 
enhanced by a factor 10 if a merging event takes place. The 
duration of the period of high efficiency star formation, tj, 
scales with the dynamical time o f the galaxy and is another 
param eter of the model. (See also lBalland. Devriendt fc Silkl 
I2OO2I) 

2.2 The space density of starbursting galaxies at 
high redshift 

Fig. shows the space density of galaxies as a function of 
star formation rate at a series of different redshifts. Overall, 
the typical star formation rates increase towards higher red- 
shift. This is because the dynamical times of high-redshift 
galaxies are shorter, implying higher star formation rates, 
and because galaxies are typically more gas-rich at high z. 
However, the space density of galaxies with star formation 
rates above 50 - 100 M0yr~^ drops at redshifts greater than 
2-3. Such systems only occur in the most massive dark mat- 
ter haloes, which are rare at early epochs. Objects with star 
formation rates of a few hundred to a thousand M0yr~^ are 
very rare at all epochs. 

To demonstrate that this is a genuine feature of all hi- 
erarchical galaxy formation models and not an artefact of 
the particular star formation and feedback prescriptions we 
have chosen, we have carried out the following simple exer- 




Figure 2. The co-moving density of galaxies that produce stars 
at rates 10^ and 10'^ MQyr~^ as obtained directly from the Press— 
Schechter mass function, t, is the star formation timescale, which 
we have related to the Hubble time, as shown in the figure. 



Figure 3. Upper diagram: the SED of a sta rburst in a giant 
molecular cloud calculated with the model of lEfstathiou et all 
(200(f). The dashed line is a grey-body spectrum fitted to the 
section of the SED with A > 100 ^m. The parameters are Ts = 34 
K, and f3s = 1.9. Lower diagram: the SED of a cirrus galaxy 
calculated with the model of ERR. The parameters are Ts = 22 
K, and f3s = 1.9. 



cise. We assign a fiducial star formation rate to each dark 
matter halo in the Universe by assuming that at any given 
epoch, one tenth of the available baryons in the halo are 
forming stars over a timescale t*. The star formation rate 
per halo may then be written as M = (n<,/r 2o)(A/fe /^«), 
wher e fJ, ris/lO, D.b = 0.02ft"^ llBania. Rood fc Balseil 
|2002^ and fio = 0.3 {Qb and Oo are the baryonic and to- 
tal co-moving mass densities in the Universe and Mh is the 
mass of the halo harbouring the galaxy). We consider two 
cases: quiescent star formation with t, — in, and bursty 
star formation with = in/ 100 and calculate halo abun- 
dances using the standard Press-Schechter formula. Fig. |5| 
shows the redshift evolution of the space density of objects 
with different star formation rates. As can be seen, the space 
density of massive haloes with ~ 10'^ Moyr"^ is small and 
dro ps very rapid l y tow ards larger redshifts. 

IScott et al.l J2Qq3) derived a density <1>(M > 1000 
M0yr-\z > 1.5) ~ lO"'^ Mpc"^ for SCUBA sources. Figs. 
and 1^ demonstrate that in 'realistic' models such space 
densities are only achieved for objects with much lower star 
formation rates of ~ 10^ M0yr~^. The common solution 
of making the star formation more 'bursty' does not help 
us match these high space densities unless star formation 
timescales are considerably shorter than we have assumed. 
This argues in favour of lower star for mation rates of the 
SCUBA sources than those inferred bv IScott et alJ (|22q2) 
unless the star formation efficiency in the SCUBA sources 
is extremely high. 



3 SPECTRAL ENERGY DISTRIBUTIONS 

3.1 Grey-body models of the Far-Infrared SED 

For a given illuminating radiation field, optical depth and 
age, the spectral energy distribution of a galaxy in the far- 
infrared (FIR) and sub-mm part of the spectrum is deter- 
mined by the dust temperature and grain size. This part of 
the S E D is relative l y featu reless (Fioc & Rocca-Volmerange 
19971: ISilva et alj Il998l: IPevricndt. Guiderdoni fc Sadat 
199£ ; l^fe^^hio^^!owan- Robinson fc Siebenmorgenn200oll 7 
because in clouds that are optically thick to optical/UV ra- 
diation, the dust emission can be qualitatively represented 
by a single large grain model whose extinction efficiency 
changes with the inverse of the grain size. The presence of 
grains of smaller size results in a distortion of the grey-body 
SED of the single grain model. We write the grey-body SED 
as Lviy) oc ei,Bi,{T,iy), where is the spectral luminosity 
of the galaxy, B,y is the Planck function, and £v{o(. i^^) is the 
dust emissivity function. In the simple, single-component 
dust model, (3 — 1. The presence of smaller grains results 
in /3 > 1. Provided that the grain size distribution is not 
changed through creation and destruction of different-sized 
grains, /3 depends only weakly on temperature. In this paper, 
P is taken as a free parameter of the model and the effec- 
tive temperature of the dust, T, is the second adjustable 
parameter. 

Fig. Bl shows g rey-body fits to the SEDs of 
lEfstathiou et all i2000h and ERR in the wavelength range 
A > 100 fj,m. The fit to the FIR emission is indeed very 
good. We therefore take a simple grey-body parameteriza- 
tion for the SED to predict the counts from our galaxy for- 
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mation model. We have compared counts calculated using 
SEDs from full radiative transfer models with those calcu- 
lated from grey-body fits. The results are in agreement to 
within ~ 1 per cent. 

The next step is to determine the normalization of the 
SED. We adopt the common practice of relating the bolo- 
metric luminosity of a galaxy to its instantaneous star for- 
mation rate through a linear relationship, ignoring the cu- 
mulative effect of old stars, and thus the star formation his- 
tory of the source. We do this in order to keep the number of 
dust parameters small. The form of the assumed initial mass 
function (IMF), the evolution of the metallicity and the dust 
mixture and the timescale of star formation all contribute to 
the uncertainty in the proportionality constant of the model, 



Lboi = CM. 



(2) 



In the literature C ^ 10^ to 10'^° Lcri M^^yr for a range 
of IMFs (see review bv lBlain et"al]l2002ll . Here, we adopt a 
value of 5x 10^ and note that it is uncertain by a factor of two 
in either direction. We further assume that the galaxies are 
optically thick at visible and UV wavelengths and that most 
starlight is absorbed and re- radiated by dust, i.e. Lgr ~ I/boi. 

It is customary to allow the high frequency limit of the 
slope of the SED to be another parameter of a grey-body 
SED model. The effect of this parameter in our calculation 
is very small since we are interested in the flux far beyond 
the FIR peak. This parameter only affects the normalization 
of the SED and its effect is a change in of < 10 per cent 
or a change of ~ 1 K in the temperature. We shall therefore 
omit it in the interest of keeping the number of parameters 
small. 



3.2 The sub-mm counts within hierarchical 
models with a hot starburst-like SED 

In ^'2.2\ we saw that, in hierarchical galaxy formation mod- 
els, the predicted space density of galaxies with star forma- 
tion rates > 100 M0yr~^ is small. This is the main reason 
for the difficulty encountered in fitting the sub-mm source 
counts when assuming a hot SED. Higher dust temperatures 
move the peak of the SED towards shorter wavelengths. 
Thus, for a given bolometric luminosity of the source, the 
hotter the SED, the smaller the proportion of the galaxy's 
output in the sub-mm. 

To demonstrate this more clearly, the curves in Fig. 
13 show the sub-millimetre counts predicted for the semi- 
analytic models, assuming that all galaxies have hot SEDs. 
Results are shown for different values of the dust emissivity 
parameter, (3. Even though the temperature is at the lower 
end of the frequently assumed range (40 - 60 K), the pre- 
dicted counts fall significantly short of the observed values. 
Increasing either the temperature or the emissivity index 
makes matters worse. 

We saw in H2.2l that the abundance of sources with mod- 
erate star formation rates of the order M ~ 50 - 100 MQyr~^ 
is sufficient to match the number of objects detected in sub- 
mm surveys. For these moderate star formation rates, it is 
unnecessary to assume a hot SED. In the next section, we ex- 
plore the effect of varying the dust temperature more freely, 
while in H5.ll we will discuss the implications for the relation 
of sub-mm sources to Lyman-break galaxies. 



10* 



in 
A 




♦ Scott et al. 3001 

O Eales et al. 2000 

□ Barger et al, 1999 

X Blaln et al. 1999 

O Hughes el al. 199B 

A Small et al. 1997 



mJv 



Figure 4. Integral number counts at A = 850 fim. Hierarchical 
models under-predict the counts at the bright end if a hot SED 
is chosen for a starburst. C, = 5.0 X 10® LQMg'^yr, tf, = 3 X 10® 
yr. In the shaded area around the solid line, ^ is varied between 
2.0 X 10® and 10^" L^M;: 



3.3 



Temperature as a free parameter in models 
with diffuse extended dust 



In the literature, a range of temperatures has been used 
when fitting sub-mm counts, but most authors have re- 
frained from allowing T to vary outside the range 40 - 60 K, 
arguing that temperatures much below this are unrealistic 
in a starburst galaxy, where star formation takes place in 
optically thick giant molecular clouds (GMCs). In standard 
radiative transfer models, the dust temperature of a GMC 
is determined by its physical extent, which is in turn limited 
by the time before it is dispersed by supernovae. The bolo- 
metric luminosity of a galaxy then scales with the number 
of molecular clouds contained within it. 

For galaxies with moderate star formation rates, it is 
plausible that most of the IR emission is due to a different 
geometrical configuration of the star/dust distribution, more 
similar to that responsible for the cirrus emission in spiral 
galaxies (ERR). If the starlight is absorbed by diffuse ex- 
tended dust and re-radiated in the FIR, the intensity of the 
UV radiation heating the dust is not set by the properties 
of GMCs and much lower dust temperatures are expected. 
The peak in the FIR emission then moves much closer to 850 
/xm and the resulting bolometric correction and the inferred 
star formation rate decrease strongly. Fig. |K| shows how the 
inferred star formation rate for a source with a ffux of 1 mjy 
at 850 fim varies as a function of redshift for a range of SEDs 
corresponding to different dust temperatures. 

3.4 Further constraints on the dust models 

The 850-/im source counts are not able to constrain all the 
free parameters needed to predict the sub-mm emission from 
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Figure 5. The star formation rate of a source with a flux of 1 
mjy at 850 fim for a range of hot and cold SEDs as a function of 
redshift. 

the galaxies in our model. We have thus chosen plausible 
values for two of our parameters, — 5.0 x 10^ LeMg'^yr 
and tb = Sx 10** yr. 

One way to break the degeneracy between star forma- 
tion rate and temperature is to consider the inferred size of 
the emitting region (Rowan- Robinson 200lj). The predicted 
radius of the emitting region depends on its bolometric lu- 
minosity and its surface brightness as 

47vr^ (3) 

J fir 

Surface brightness and temperature are linked because of 
the equilibrium between absorbed UV/optical emission and 
emitted IR emission, 1^^ oc /heat- Higher UV/optical sur- 
face brightness, /heat, will result in higher dust temperatures 
which in turn will result in a higher IR surface brightness. 

In order to normalize the relation between the intensity 
of the radiation field and the emitted bolometric surface 
brightness, we have to make assumptions about the relative 
distribution of dust and stars. Observed sub-mm sources are 
generally very faint in the optical and UV and are therefore 
heavily obscured at these wavelengths. We will consider two 
limiting cases: a dust screen and an extended distribution 
of stars and dust. For a dust screen, the situation is basi- 
cally the same as in the optically thin case. The dust layer 
located within about one optical depth from the stars will 
be heated to approximately the same equilibrium tempera- 
ture as in the optically thin case. The surface brightness of 
starlight needed to reach a certain temperature is fixed and 
can be used to infer the bolometric FIR surface brightness 
and hence the size of the IR-emitting region. 

The dust properties (composition, grain size distribu- 
tion) will fix the dust emissivity index f3. The predicted 
surface brightness/size will increase/decrease with increas- 
ing star formation rate. The relations between tempera- 



Figure 6. Integral number counts at 850 lira. The sources of 
the data are given in the legend. The error bars in each data set 
are correlated. ^ = 5.0 X 10^ L0M~^yr, tj, = 3 X 10* yr unless 
otherwise stated. 

ture, bolometric luminosity, bolometric surface brightness 
and sub-mm flux can then be used to predict the sizes of 
sources of given sub-mm flux in our models, or, more gener- 
ally, to infer a star formation rate from an observed sub-mm 
flux and size. We will discuss this in more detail in ^4.'2\ 



4 RESULTS 

4.1 Fitting the 850-/im counts 

In this section, we use the star formation rates given by the 
hierarchical galaxy formation model, described in ^21 to cal- 
culate the 850-//m counts. For simplicity, we have assumed 
that the sub-mm sources are all star-forming objects char- 
acterized by a single temperature SED. Of course, in real- 
ity, sources of a range of temperatures, most likely also be- 
longing to different populations of sources (quiescently star- 
forming galaxies, starbursts, AGN, etc.), will contribute to 
the 850-/im counts. 

There is a certain amount of degeneracy between the 
temperature and the emissivity index. As illustrated in Fig. 
|S] the temperature required to obtain a reasonable fit to 
the data changes from 23 K to 21 K as we vary /3 in the 
range 1.65 to 1.95. In physical terms, the value of l3 is set 
by the precise nature and composition of the dust mixture, 
its optical depth, and how these vary as a function of the 
age of the stars that are re sponsible for the heating. Realistic 
radiative transfer models iSilva et alll99d : lEfstathiou et all 
^OOft ERR) yield l3 values in the range 1.5 - 2.0. 

The parameter ^ parameterizes the effect of the IMF on 
the energy output by a given mass of stars that is formed. 
A change in moves the count distribution horizontally. A 
change of the star formation time scale, tb, alters both the 
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luminosity of the sources and the length of time they are 
'visible', so its effect on the count distribution is more com- 
plex. The bright counts are dominated by the most strongly 
star-forming objects. Because the star formation rate func- 
tion is very steep at the bright end (see Fig.0, an increase 
in the luminosity of the sources tends to over-compensate 
for the drop in density. The faint end of the star formation 
rate function is shallow, so at fainter flux levels the normal- 
ization of the counts is determined by the space density of 
the sources. The cross-over of the two regimes takes place 
at ~ 0.1 mjy. Over the range of observed source fluxes, a 
lower value of will thus reduce the slope of the counts. 

The effect of a scatter in the temperatures of the sub- 
mm emitting regions is to reduce the slope of the flux dis- 
tribution, as lower temperature sources will make a large 
contr ibution to the counts at the bright end iFardal et alJ 
I2OOII). 



4.2 Sizes of the Source Galaxies 

Because of its large beam size, SCUBA is unable to constrain 
the sizes of the detected sub-mm sources. The only informa- 
tion we have so far on the sizes of these sources comes from 
observations with the Plateau de Bure interferometer which 
appear to have marginally resolved the SCUBA galaxy Lock- 
man850.1 (Lutz ct al. 2001). As we will discuss below, the 
size of this galaxy fits in with the assumption of low dust 
temperatures. 

Let us define the SED of a source as Lv — fvLtoi, where 
is a normalized function. We can then use the fiux redshift 
relation to relate flux and luminosity. 



(4) 



Here, Dl is the luminosity distance. Equation |^ then gives 
an expression for the angular size of the object in terms of 
its bolometric surface brightness Jfir, 



QS, 



-1/2 



-1/2 



(5) 



As discussed in H8.4I for a given grain composition, the 
balance between absorbed and emitted energy establishes a 
tight relation between the intensity of the UV/optical radi- 
ation and the temperature (distribution) of the dust. In the 
case of an extended, optically thin dust distribution, this in- 
tensity will be given by the average surface brightness of the 
starlight (as is the case for cirrus emission from spiral galax- 
ies). As discussed, the resulting FIR SED can be well fitted 
by a grey-body spectrum with an effective temperature and 

T^+^ OC /heat. 

We have used the cirrus model of ERR to calculate the 

SED and the relation between temperature and /heat ■ If the 

dust acts as a screen /hr ~ (1 — e""*^ )/heat ^ /heat, while for 

an extended distribution of dust and stars /gr ~ ^v/hoat- 

Fig. 13 shows the resulting angular sizes as a function of 

redshift for sources with S^^ — 1 mJy, assuming an extended 

dust distribution with Ay = 10. Results are shown for a 

range of different dust temperatures. The value of Av we 

assume is slightly larger than that adopted by ERR to fit the 

SEDs of non-starburst SCUBA sources. However, in order to 

have optical data, the sources they consider are necessarily 

1/2 

less obscured. Note the obvious scaling oc S,^g . For a dust 



<s 




Figure 7. The angular size of a source with a flux of 1 mJy at 
850 fim for the fitting models and for the limiting cases in the 
commonly assumed range of 40 K < T < 60 K and 1 < /3 < 2. 
The data point is taken from lLutz et al ]|200i) 



screen geometry, the predicted sizes would be larger by a 
factor three. 

The angular size of the sources is almost constant at 
low redshifts. Once the k-correction no longer compensates 
for the redshift factor in equation |S] the angular sizes in- 
crease with redshift. For low dust temperatures, the sizes 
are around l(S'85o/lmJy)^''^ arcsec for redshifts less than 
~ 2 and then increase by a factor 2 in the redshift interval 
2 < z < 6. We also show two models with temperature and 
emissivity appropriate for hot starburst SEDs. At a given 
redshift, the angular size of the source decreases approxi- 
mately linearly with the assumed temperature. The cool, 
extended models predict significantly larger sizes than the 
hot SED models. A galaxy with {T, 13) = (22 K, 1.8) would 
be two to three times as large as a hot source in the redshift 
range 2 < z < 6. 

The weak redshift dependence of the angular sizes 
makes it possible to use observed sizes to constrain the dust 
temperature in the screen model and a combination of tem- 
perature and optical depth for the extended model. 

The cross in Fig. [7| shows the size of SCUBA galaxy 
Lockman850.1 assuming that it is marginally resolve d and 
scaled to an 850- ^im fiux of 1 mJy jLutz et alJl200lD . This 
observation is in reasonably good agreement with our pre- 
diction for emission from cold extended dust. Note, however, 
that LE850. 6, resolved in the r adio seems to have a slightly 
smaller size Jivison et alJl2002ll . 

The corresponding proper sizes of the sources depend 
on the cosmology. With the same parameters as in Fig. |S| 
our model predicts sizes > 5 kpc (Fig.|HJ. 

Changes in and it require changes in the assumed 
dust temperature and emissivity index in order to maintain 
a good fit to the counts. An increase in means that the 
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Figure 8. The size of a source with a flux of 1 mjy at 850 
for the fitting models and for the limiting cases in the commonly 
assumed range of 40 K < T < 60 K and 1 < /3 < 2. 



Figure 9. Redshift distribution of the sources for the model with 
C = 5 X 10^ LqMq Vr and = 10^ Gyr and T = 22 K, at S^^ = 1 
mJy (solid), 0.3 mJy (dotted) and 10 mJy (short-dashed). The 
other two lines are both at 1 mJy and = 10^ Gyr, = 5 X 10^ 
LqUq^jv, T = 17 K (long-dashed) and C = 8.9 X lO'' LoMgVr, 
T = 21 K (dot-dashed). All have /3 = 1.8 



dust temperatures can be higher, since a smaller fraction of 
the stellar output is needed in the sub-mm. Thus the sources 
can be more compact, with larger surface brightnesses and 
smaller radii. Increasing the star formation timescale has 
the opposite effect. The bolometric source luminosities drop 
and cooler, more extended sources are required to match the 
observed sub-mm fluxes. 



4.3 Redshift Distributions 

Fig. |S| shows the redshift distribution of our sources down 
to various flux limits for a range of different parameters. 
The typical redshifts are in the range 2 < z < 4. This is 
very similar to previous modelling which has assumed hot 
starburst-like SEDs. 

The solid curve shows the redshift distribution down 
to a flux level of 1 mJy with the parameter set (T, /3) = 
(22 K, 1.8). The median redshift of the sources is 2.6. The 
mean of the redshift distribution is 2.8 and its variance 1.4. 
The distribution has a high redshift tail that extends to 
z « 6. The distribution of the fainter sources is similar, but 
extends to larger redshifts (dotted) while the bright sources 
peak at slightly lower redshifts (short-dashed). The curves 
change little if we change the starburst time scale and the 
assumed IMF (long-dashed and dot-dashed, respectively) 

Because the sub-mm counts are strongly biased towards 
high-redshift sources, counts at shorter wavelengths place 
rather little constr aint on models for t he sub-mm source 
population (see e.g. iGranato et al]|200(J) . 



5 DISCUSSION 

5.1 The nature of sub-mm sources and the 
relation to Lyman— break galaxies 

As discussed in the previous section, low dust temperatures 
imply a star/dust distribution that is spatially extended over 
> 5 kpc. This is very different from typical violent starbursts 
at the centre of galaxies which are spatially compact with 
typical sizes < 100 pc. In our model, the sub-mm sources 
are more extended than Lyman-break galaxies which have 
typical sizes of 0.6 arcsec liCalzetti fc GiavaUscollioOlh . 

The UV continuum slope of Lyman-break galaxies in- 
dicates that they are typic ally optically thick in the opti- 
cal/UV JPettini et alj|200 j) and emit the bulk of their Ught 
in the IR. Typical correction factors to the bolometric lu- 
minosity rang e from 3 to 15 depending on the assumed red- 
dening curve IS teidel et al ]|W98). Furthermore, the spectra 
of Lyman-break galaxies exhibit many of the characteristics 
of nearby starbursts. It was thus expected that they would 
dominate the sub-mm source population and the lack of 
sub-mm detections of most Lyman-broak galaxies came as a 
surpr ise ( Hughes et al. 1998; Sniail ct a l. 1998; Bargor et aO 
ll999l:lmTerariri999l:IChapman et allboOOl) . 

The larger spatial extent and smaller surface bright- 
ness predicted by our cold dust models may naturally ex- 
plain the lack of overlap of the sub-mm sources and Lyman- 
break galaxies. Optical/UV selected samples of high-redshift 
galaxies will be strongly biased to high surface brightness 
objects. These will have high dust temperatures and, at the 
typical redshifts of these sources, the 850-/im flux will come 
from far down the Raleigh- Jeans tail of the spectrum. Con- 
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versely, a larger fraction of tlie total emitted light of low 
surface brightness objects will be detected in the sub-mm. 
In this way, the UV to sub-mm flux ratio can differ by a 
factor 30 - 100 for objects with the same optical depth and 
star formation rate (see Fig. 01. A larger optical depth in 
sub-mm detec ted sources would obv iously further increase 
this ratio (see IShu. Mao fc M0II2OOII for a model where the 
bright sources are more strongly obscured). 

5.2 Contribution from AGN 

Active Galactic Nuclei (AGN) make a contribution to the 
FIR and sub-mm luminosities of Luminous IR Galaxies 
by heating the d ust in torus enshrouding an AGN (e.g. 
lRowan-Robinsonl lT995l . The debate on the fraction of the 
infrared luminosity due to AGN has not yet been set- 
tled, though recent surveys with Chandra seem to sug- 
gest that the co ntribution of AGN to SCUBA sources 
is not signiflcant l| Alexander et all 1200 J) . Observations of 
galaxies in the local universe suggest 70 - 80 per cent 
of ULIRGs (lO^Lffl < Lfir < IO^^Lq) are starburst- 
dominated iGenzel et alJ Il998l) . Only when the luminos- 
ity reaches that of a Hyper-luminous IR Galaxy (HLIRG, 
Lfir > lO^^ Lpi), are AG N thought to dominate the far- 
IR emission llFarrah et alj E)02'l. For almost all the sources 
m our model Lir/Lq < 10^^ ^ The average contribu- 
tion from AGN to these galax ies is of order 1 per cent 
JVeilleux. Kim fc Sander3ll999|) . 

5.3 Implications for source counts in the mid IR 

The forthcoming Spectral and Photometric Image Receiver 
(SPIRE; Griflin, Swinvard & Vigroux 2001) on European 
Space Agency's Herschel Space Observatory will image the 
sky at 250, 350, and 500 /^m. The high redshifts of our source 
population means that the first of these channels already 
looks at emission from beyond the FIR peak of our SEDs 
{z ~ 1.2). Only the brightest low-redshift SCUBA sources 
would be detectable in a 250 /im SPIRE survey. At 350 
Hm, SCUBA fluxes are almost unchanged for sources at the 
median redshift. The 500 /xm channel will see most objects 
with Ss50 5 mjy over the whole redshift range. The point 
source sensitivity of the channels is adequate to measure 
the slope of the SED, and hence the temperature, of at least 
the brightest of the sources, directly. If sources exhibit some 
scatter in temperature (as is likely), all the above surveys 
will be biased toward the sources with higher temperatures 
than we have assumed in this paper. 



6 CONCLUSIONS 

The sub-mm source counts predicted by hierarchical models 
of galaxies formation have been studied for a wide range of 
SEDs. 

We confirm the results of previous studies that, with 
plausible assumptions about star formation timescales, star 
formation efficiencies and IMF, hierarchical models for 
galaxy formation under-predict the number of bright sub- 
mm sources by a large factor if hot (40 - 60 K) starburst-like 
SEDs are assumed. 

The sub-mm source counts alone do not constrain the 



dust temperatures and SEDs of the sources. Dust temper- 
atures and inferred star formation rates are highly degen- 
erate. Lower dust temperatures require substantially lower 
star formation rates to produce the same 850-/im flux. 

For dust temperatures in the range 20 - 25 K, the ob- 
served sub-mm counts agree very well with those predicted 
by hierarchical galaxy formation models. These low tem- 
peratures imply typical radii that correspond to angular 
sizes of l(S'85o/lmJy)^''^ arcsec if we assume an extended 
distribution of dust and stars with Ay = 10. This is a 
factor two to three larger than predicted for starburst-like 
SEDs. Observational limits on the size of sub-mm sources 
are scarce. Marginally resolved observations of the SCUBA 
galaxy Lockman850.1 with the Plateau de Bure interferom- 
eter, and follow up near infrared imaging, may indicate that 
the sub-mm sources are as large as predicted by our cold 
extended dust models. The typical redshift of the sources 
should then lie in the range 2 < 2: < 4, somewhat, but 
not signiflcantly, smaller than in models which assume a 
starburst-like SED. These low temperatures substantially re- 
duce the predicted mid-IR emission of the sub-mm sources. 
If our models are correct, SCUBA-selected sources may con- 
tribute less than expected to the counts in upcoming mid- 
IR surveys, which will be dominated by sources with hotter 
dust. 

The low temperatures required to reconcile the sub-mm 
counts with the predictions of hierarchical galaxy formation 
models are not plausible for compact star-bursting galaxies, 
but they are expected if the emission comes from a more 
extended distribution of stars and dust. If most sub-mm 
galaxies have indeed extended distributions of cold dust with 
large sub-mm to UV/optical ratios this would naturally ex- 
plain the small overlap with the population of Lyman-break 
galaxies, which is biased towards high surface brightness re- 
gions in actively star-forming galaxies. 
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